Synchronization and bistability of qubit coupled to a driven dissipative oscillator 
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We study numerically the behavior of qubit coupled to a quantum dissipative driven oscillator 
(resonator). Above a critical coupling strength the qubit rotations become synchronized with the 
oscillator phase. In the synchronized regime, at certain parameters, the qubit exhibits tunneling 
between two orientations with a macroscopic change of number of photons in the resonator. The 
life times in these metastable states can be enormously large. The synchronization leads to a 
drastic change of qubit radiation spectrum with appearance of narrow lines corresponding to recently 
observed single artificial-atom lasing [O. Astafiev et al. Nature 449, 588 (2007)]. 
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In physics there are not so many simple quantum prob- 
lems which are exactly solvable [1]. Two of them are 
monochromatically driven two-level atom (spin-half or 
qubit) [2] and quantum oscillator (unitary or dissipative) 
[3, 4]. One atom weakly coupled to a field in a resonator 
is known as the Jaynes-Cummings model which is also in- 
tegrable [5, 6, 7]. At strong coupling the dynamics may 
become nontrivial with the emergence of classical [8, 9] 
and quantum chaos [10] but in this case one should have 
many atoms which may absorb many photons. For one 
atom even with a strong coupling to a quantum photonic 
field the dynamics is still relatively simple due to a total 
energy balance [11]. 

This old problem of a two-level atom coupled to 
photons regained recently a significant interest due to 
appearance of long living superconducting qubits [12] 
which can be strongly coupled to a microwave resonator 
[13, 14, 15]. There are also other possibilities of supercon- 
ducting qubit coupling to a quantum oscillator [16, 17]. 
The oscillator can be realized as a tank circuit tuned to 
the Rabi frequency [16] or as a current-biased dc SQUID 
[17] allowing efficient energy exchange with a qubit. Pos- 
sibilities of qubit coupling to a cooled nanomechanical 
resonator are actively discussed [18, 19] and coupling be- 
tween micro-mechanical cantilever and atomic spin found 
impressive experimental implementations (see [20] and 
Refs. therein). However, the most intriguing way seems 
to be the coupling with a microwave resonator where 
a lasing has been realized recently with 6-30 photons 
pumped into the resonator [15]. Contrary to recent in- 
teresting theoretical studies [20, 21, 22] where pumping 
is applied to a spin or qubit we concentrate here on the 
case where a monochromatic pumping is applied to a dis- 
sipative oscillator (resonator). Such a dissipative oscilla- 
tor can be also viewed as a semiclassical detector which 
performs monitoring of a qubit. This continuous type 
of measurements is now actively discussed for supercon- 
ducting qubits and other sofid-state devices [23]. The 
continuous measurement of a superconducting qubit is 
realized in [16]. 



The Hamiltonian of our model reads 

H — hujoh — hQax/2 + g7ki;o(a + a^)o'z + / cos ujt (a -f a^) 

where g is a dimensionless coupling constant, the driving 
force amplitude and frequency are / = Ti\,Jn^ and cj, 
the oscillator frequency is ojq and Ml is the qubit energy 
spacing. As in [22] we choose a qubit coupling via CTz. We 
assume that the qubit life time is enormously long and 
that its dynamics is perturbed only by the coupling with 
the driven dissipative oscillator. The dissipation rate of 
oscillator is A and we assume the quality factor to be 
Q = luq/X 100. The evolution of the whole system is 
described by the master equation for the density matrix 
p which has the standard form [4]: 

i - 11 
P = -■^[■^' P\ + Hapa'' ~ -^a^pa - -pa^a) (2) 

The numerical simulations are done by direct solution of 
time evolution of p expanded in a finite basis of oscillator 
states n, by the state diffusion method [24] and by the 
method of Quantum Trajectories (QT) [25]. We ensured 
that these methods give the same results but the major- 
ity of data is obtained with quantum trajectories which 
we found to be more suitable for massive simulations. In 
addition the QT have an advantage of providing a pic- 
torial illustration of individual experimental runs. The 
numerical details are the same as in [26, 27] and we use 
here up to n = 70 oscillator states which give good nu- 
merical convergence. Our results show that a coupling of 
two simple integrable models gives a nontrivial interest- 
ing behavior. 

A typical example of QT is shown in Fig. 1. It 
shows two main properties of the evolution: the os- 
cillator spends a very long time at some average level 
(n) = n_ and then jumps to another significantly differ- 
ent value n+. At the same time the polarization vector 
of qubit ^ defined as ^ = Tr{p(j) also changes its orien- 
tation direction with a clear change of sign of ^j, from 
S^x > Q to £,x < 0. The time averaged values of £_y^z 
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FIG. 1: (color online) Bistability of qubit coupled to a driven 
oscillator with jumps between two metastable states. Top 
panel shows average oscillator level number (n) as a func- 
tion of time t at stroboscopic integer values ujt/2n; mid- 
dle panel shows the qubit polarization vector components 
(blue/black) and (green/gray) at the same moments of 
time; the bottom panel shows the degree of qubit polarization 
^. Here the system parameters are X/u>o = 0.02, uu/uo — 1.01, 
n/wo = 1.2, / = hXy/n^, Up = 20 and g = 0.04. 



are zero but when they are taken at stroboscopic integer 
moments Lot/2'K they also show transitions between two 
metastable states. The transition time is approximately 
tm ^ 1/A being rather small compare to the life time 
in a metastable state. Inside such a state the degree of 
qubit polarization ^ = |^| is very close to unity showing 
that the qubit remains mainly in a pure state. The drops 
of ^ appear only during transitions between metastable 
states. Special checks show that an inversion of by 
an additional pulse (e.g. from > to < 0) pro- 
duces a transition of oscillator to a corresponding state 
(from n_ to n+) after time tm ~ 1/A. Thus we have here 
an interesting situation when a quantum flip of qubit 
produces a marcoscopic change of a state of detector (os- 
cillator) which is continuously coupled to a qubit (we 
checked that even larger variation n± ^ Hp is possible by 
taking Up = 40). In addition to that inside a metastable 
state the coupling induces a synchronization of qubit ro- 
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FIG. 2: (color online) Top panels: the Poincare section 
taken at integer values of ujt/2TT for oscillator with x — 
{{a + a'^)/V2), {{a- a})/y/2i) (left) and for qubit polar- 
ization with polarization angles {9,(f)) defined in text (right). 
Middle panels: the same quantities shown at irrational mo- 
ments of uot/2-K. Bottom panels; the qubit polarization phase 
cj> vs. oscillator phase {p/x = —taiitf) at time moments 
as in middle panels for g = 0.04 (left) and g — 0.004 (right). 
Other parameters and the time interval are as in Fig. 1. The 
color of points is blue/black for > and red/gray for 

^0= < 0. 



tation phase with the oscillator phase which in its turn is 
fixed by the phase of driving field. The synchronization is 
a universal phenomenon for classical dissipative systems 
[28]. It is known that it also exists for dissipative quan- 
tum systems at small effective values of h [27]. However, 
here we have a new unusual case of qubit synchronization 
when a semiclassical system produces synchronization of 
a pure quantum two-level system. 

The phenomenon of qubit synchronization is illus- 
trated in a more clear way in Fig. 2. The top panels 
taken at integer values ujt/2TT show the existence of two 
fixed points in the phase space of oscillator (left) and 
qubit (right) coupled by quantum tunneling (the angles 
are determined as = ^cos6,^y = ^sin^sin^, — 
£, sin 9 cos (j>) . A certain scattering of points in a spot of 
finite size should be attributed to quantum fluctuations. 
But the fact that on enormously long time (Fig. 1) the 
spot size remains finite clearly implies that the oscillator 
phase if is locked with the driving phase ut inducing the 
qubit synchronization with (p and ujt. The plot at t val- 
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FIG. 3: (color online) Right panel: dependence of average 
qubit polarization components and (full and dashed 
curves) on g, averaging is done over stroboscopic times (see 
Fig.l) in the interval 100 < ijjt/2-K < 2 x 10''; color is fixed by 
the sign of averaged over 10 periods (red/gray for < 0, 
blue/black for > 0; this choice fixes also the color on right 
panel). Left panel: dependence of average level of oscillator 
in two metastable states on coupling strength g, the color is 
fixed by the sign of on right panel that gives red/gray for 
large n+ and blue/black for small n_; average is done over 
the quantum state and stroboscopic times as in the left panel; 
dashed curves show theory dependence (see text)). Two QT 
are used with initial value = ±1- AH parameters are as in 
Fig. 1 except g. 




FIG. 4: (color online) Dependence of average level n± of os- 
cillator in two metastable states on the driving frequency uj 
(average and color choice are the same as in right panel of 
Fig. 3); coupling is g = 0.04 and g = 0.08 (dashed and full 
curves). Inset shows the variation of position of maximum at 
UJ = u>± with coupling strength g, Auj± — uj± — ujo- Other 
parameters are as in Fig. 1. 



ues incommensurate with 27r/w (middle panels) shows 
that in time the oscillator performs circle rotations in 
{p, x) plane with frequency lo while qubit polarization 
rotates around a:-axis with the same frequency. Quan- 
tum tunneling gives transitions between two metastable 
states. The synchronization of qubit phase with oscil- 
lator phase ip is clearly seen in bottom left panel where 
points form two lines corresponding to two metastable 
states. This synchronization disappears below a certain 
critical coupling gc where the points become scattered 
over the whole plane (panel bottom right). It is clear that 
quantum fluctuations destroy synchronization for g < gc- 
Our data give gc — 0.008 for parameters of Fig. 1. 

The variation of bistable states with coupling strength 
g is shown in Fig. 3. The difference between ri+ and 
n_ grows with g. It is striking that n_ may become 
close to zero. The direction of qubit polarization also 
changes in a smooth but nontrivial way. It is also im- 
portant to note that according to our data the disper- 
sion of oscillator wave function in metastable states is 
compatible with the dispersion of a coherent state with 
n±. This corresponds to a wave packet collapse induced 
by dissipation (see [26, 27] and Refs. therein). The 
dependence of n± on driving frequency uj is shown in 
Fig. 4. A symmetric double peak structure is evident: 
for oj > wq the metastable state with < has maxi- 
mal n value while for u; < loq the state with maximal n 
has > (note color interchange). The peak width is 
approximately equal to the dissipation rate A. With the 
increase of g their form becomes asymmetric indicating 
importance of nonlinear effects. The splitting of peaks 
grows approximately linearly with g (inset at Fig. 4) and 
reminds the vacuum Rabi splitting effect [6]. The shift 



Auj± explains two states n± of driven oscillator well de- 
scribed by n± = npX'^/{4:{uj — ivq — Auj±)^ + A^)] (see 
dashed curves in Fig. 3 left traced with numerical values 
of Auj± from Fig. 4 inset). To estimate Auj± we note that 
the frequency of effective Rabi oscillations between quasi- 
degenerate levels is ilu ~ gi^Q\^n± + 1 [5, 7] that gives 
Au!± ~ (Klu/dn sa ±gujQ/2^n± -f 1 in a good agreement 
with data of Fig. 4 for moderate g. 

The properties of two metastable states are analyzed 
in Fig. 5. The number of transitions Nf between these 
states has a pronounced peak at f2 « l.lwo that approx- 
imately corresponds to a resonance condition Vl — uj k, 
2guJo. For fl < l.OSwo there is an abrupt drop of Nf and 
bistability becomes irregular disappearing for certain 
but the synchronization still remains. Quite interestingly, 
the data show that for fl > 1.1 the life times of each 
state are rather different and enormously large, generally 
T_ > r+ » ojq/X. 

The spectrum of qubit radiation S{i') in presence of 
phase noise in (jj is shown in Fig. 6. It confirms the main 
features discussed above: for Q/luq — 1.2 the growth of 
driving power Up induces the synchronization of qubit 
with radiation suppression at qubit frequency CI = 1.2^0 
and appearance of narrow line with lasing at v — lu. For 
il = ujo the radiation spectrum S{i') at rip < 1 has two 
broad peaks at v = uj± gtoo corresponding to the vacuum 
Rabi splitting [6] (narrow line from driving source at v = 
UJ is also visible in this case). At strong driving np > 1 
the synchronization takes place with appearance of one 
lasing line at v = uj. For both values of Q the transition 
to synchronization/lasing takes place at Hp > npi « 2. 
The spectrum S{v) in Fig. 6 has close similarities with 
the spectrum observed recently in a single artificial-atom 



4 




FIG. 5: (color online) Dependence of number of transitions 
Nf between metastable states on rescaled qubit frequency 
Sl/wo for parameters of Fig. 1; Nf are computed along 2 QT of 
length 10^ driving periods. Inset shows life time dependence 
on Q/ll>o for two metastable states (r+ for red/gray, t_ for 
blue/black, t± are given in number of driving periods; color 
choice is as in Figs. 2,3). 
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FIG. 6; (color online) Spectral density S{v) of qubit radiation 
^z{t) as function of driving power rip in presence of phase 
noise in </!> with diffusion rate rj = 0.004a;o. Left: Q/ujo ~ 
1.2; right: i}/ujo — 1. Other parameters are as in Fig. 1. 
Color/grayness shows S{i') in logarithmic scale (white/black 
for maximal/zero), u is given in units of uo- 

lasing [15] which appears at a similar threshold Upi w 1. 
A shift related to splitting w± ~ gwo/Zy^ is also seen 
experimentally. However, the vacuum Rabi splitting is 
not visible in Fig. 3c of [15]. Exact comparison requires 
much more extended numerical simulations since in [15] 
Q w 10-* while we have Q ~ 100. 

In conclusion our numerical studies show a nontrivial 
behavior of a rather simple model given by Eqs. (1,2). 
It is characterized by bistability and synchronization of 
qubit induced by its coupling to a quantum driven dis- 
sipative oscillator. As for the vacuum Rabi splitting [6] 
it is important that the oscillator is quantum since the 
effect is absent for a classical dissipative oscillator with 
commuting a,a^ (1). A better analytical description and 
understanding of the behavior discussed requires further 
studies. Especially interesting is the analysis of long life 
times T± in metastable states related to a macroscopic 
quantum tunneling. 



Preparing the paper to submission we became aware 
of the preprint [29] where a similar model is studied but 
synchronization is not discussed there. 
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